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We propose that dark matter is not yet another new particle in nature, but that it is a remnant 
of quantum gravitational effects on known fields. We arrive at this possibility in an indirect and 
surprising manner: by considering retarded, nonlocal, and Lorentzian evolution for quantum fields. 
This is inspired by recent developments in causal set theory, where such an evolution shows up as 
the continuum limit of scalar field propagation on a background causal set. Concretely, we study 
the quantum theory of amassless scalar field whose evolution is given not by the the d’Alembertian 
□ , but by an operator □ which is Lorentz invariant, reduces to □ at low energies, and defines an 
explicitly retarded evolution: (D0)(a;) only depends on rf>(y), with y is in the causal past of x. This 
modification results in the existence of a continuum of massive particles, in addition to the usual 
massless ones, in the free theory. When interactions are introduced, these massive or off-shell quanta 
can be produced by the scattering of massless particles, but once produced, they no longer interact, 
which makes them a natural candidate for dark matter. 

PACS numbers: 


I. INTRODUCTION 

The nature of dark matter is one of the most impor¬ 
tant problems in modern physics. Almost a century af¬ 
ter it was hypothesized, though, our understanding of 
it is still limited to its gravitational signature on lumi¬ 
nous matter. It is often assumed that dark matter is 
a new weakly interacting particle which is just hard to 
detect. However, so far there has been no conclusive 
direct or indirect detection in accelerators or cosmolog- 
ical/astrophysical settings. In what follows, we propose 
that dark matter is not yet another new particle in na¬ 
ture, but that it is a remnant of quantum gravitational 
effects on known fields. We arrive at this possibility in an 
indirect and surprising manner: by considering retarded, 
nonlocal, and Lorentzian evolution for quantum fields. 
Concretely, we study the consequences of replacing the 
d’Alembertian □ with an operator □ which is Lorentz 
invariant, reduces to □ at low energies, and defines a re¬ 
tarded evolution: (□<}(>)(£) only depends on with y 
is in the causal past of x. Why is this type of evolution 
interesting, what does it have to do with quantum grav¬ 
ity, and how does it lead to a proposal for the nature of 
dark matter? 

The causal set theory approach to quantum gravity 
postulates that the fundamental structure of spacetime 
is that of a locally finite and partially ordered set p]. Its 
marriage of discreteness with causal order implies that 
physics cannot remain local at all scales. This nonlo¬ 
cality manifests itself concretely, for instance, when one 
seeks to describe the wave propagation of a scalar field on 
a causal set. It has been shown in this case that coarse- 
graining the quantum gravitational degrees of freedom 
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leads to a nonlocal field theory described by an operator 
exactly of the type □ There are reasons to suspect 

that this type of nonlocality is not necessarily confined to 
the Planck scale, and that it may have nontrivial impli¬ 
cations for physics at energy scales accessible by current 
experiments (see BIB] and references therein for impli¬ 
cations of nonlocality in the context of cosmology). It is 
then only natural to wonder what a quantum field the¬ 
ory built upon □ would look like, especially that it may 
contain information about the fundamental structure of 
spacetime. 

Studying □ is also interesting from a purely field- 
theoretic perspective, since it forces us to relax one of 
the core assumptions of quantum field theory: local¬ 
ity. Most nonlocal and Lorentzian quantum field the¬ 
ories studied in the literature consider modifications of 
the type □ —> /(□). In this paper, we consider explic¬ 
itly retarded operators, which are more generic and have 
more interesting properties as a result. For instance, the 
Fourier transform of □ is generically complex, which is 
a direct consequence of retarded evolution. In fact, this 
feature is at the heart of our proposal for the nature of 
dark matter. It is also worth mentioning that quantizing 
a field theory of the type described here is non-trivial due 
to the absence of a local action principle. This presents 
a technical challenge, from which one may gain deeper 
insight into quantization schemes. 

What is the relation between a quantum field theory 
based on □ and dark matter? Upon quantizing a free 
massless scalar field 4>{x) with the classical equation of 
motion tA(j>{x) = 0, we find off-shell modes in the mode 
expansion of the quantized field operator These 

are modes which do not satisfy any dispersion relation, 
unlike in usual local quantum field theory (LQFT) where 
every Fourier mode with four-momentum p is an on-shell 
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quanta, i.e. it satisfies p ■ p = O. 1 This is equivalent 
to the statement that the quantized field operator does 
not generically satisfy the classical equation of motion: 
□^>(x) 0. Note that an off-shell mode of a massless 

scalar field has an effective mass, and can be thought of 
as a massive quanta in itself. We show that the off-shell 
modes can exist in “in” and “out” states of scattering, 
and are different from virtual particles which exist as 
intermediate states in Feynman diagrams. When consid¬ 
ering the interacting theory, we find an extremely sur¬ 
prising result: the cross-section of any scattering process 
which contains one or more off-shell particle(s) 2 in the 
“in” state is zero. That is to say, on-shell quanta can 
scatter and produce off-shell particles, but once produced, 
off-shell particles no longer interact. It is this behaviour 
that makes these off-shell particles a natural candidate 
for dark matter. The phenomenological story would be 
that dark matter particles were produced in the early 
universe in this fashion: as off-shell modes of quantum 
fields. This feature of the theory can be traced back 
to the fact that □ defines an explicitly retarded evolu¬ 
tion, which as mentioned previously, may be a remnant 
of quantum gravitational degrees of freedom. 

Our paper is organized as follows. In Section [TTJ we 
start by setting forth a series of axioms which any non¬ 
local, retarded, and Lorentzian modification of □ at high 
energies should satisfy. In Section m we argue there 
is no action principle for the theory of interest, which 
forces us to carefully study, in Section IV what quan¬ 
tization scheme should be used. There, we argue that 
canonical quantization and the Feynman path-integral 
approach do not work, and explain why the Schwinger- 
Keldysh (also known as the double path integral or in-in) 
formalism provides the appropriate framework. Sections 
|V| and |VI| describe the interacting theory, where we work 
out the modified Feynman rules, find S-matrix ampli¬ 
tudes, and compute cross-sections for various examples 
and comment on the time reversibility of the theory. Al¬ 
though a continuum, superposition of off-shell particles 
can in principle scatter into on-shell modes, we argue 
why this is unlikely to happen. Extension to massive 
scalar fields is discussed in IVIII Section IVIIII concludes 
the paper. 


II. MODIFIED D’ALEMBERTIAN: DEFINITION 


In this section we study generic spectral proper¬ 
ties of non-local and Lorentzian modifications of the 
d’Alembertian □. We focus on a class of operators □ 


1 We use a signature of-f- ++ for the Minkowski metric 

Also, pi P2 S V^PiP2- 

2 In the quantum theory, an off-shell particle is 1 -particle quantum 
state with a well-defined (non-zero) mass and momentum, i.e. a 
massive eigenstate of Hamiltonian and momentum operator. 


which defines an explicitly retarded evolution: (□</>) (x) 
depends only on <j>(y) with y in the causal past of x. 
As we will see, such operators have interesting features 
which are absent in modifications of the type /(□). We 
start by setting forth a series of axioms which a non¬ 
local, retarded, and Lorentzian modification of □ at high 
energies should satisfy: 

1. Linearity: 

\3(a<f> + bip) = aO(f> + bOtf), a, b £ C, (1) 

where 4> and ip are complex scalar fields and C de¬ 
notes the set of complex numbers. 

2. Reality: for any real scalar field <f> , □</> is also real. 
Note that reality and linearity imply for any com¬ 
plex scalar field <p that 

(□0*) = (00)*, (2) 
where * denotes complex conjugation. 

3. Poincare-invariance: evolution defined by □ is 
Poincare-invariant. Consider a scalar field (p(x) 
which transforms to (p'{x) = </>(A _1 x) under a 
Poincare transformation x -A Ax. We require □ 
to be invariant under the action of A: 

(□(//) (x) = (□</>)(A -1 x). (3) 

Taking A to be a spacetime translation A(x) = x + 
a, one finds that the eigenfunctions of □ are plane 
waves. To see this, let <f)(x ) = e lp x and define 
ip(x) = (□</>)(x). It then follows from (J3| that 

e~ lp ' a ip(x) = ip(x — a), (4) 

where we have used the linearity condition. Solu¬ 
tions to the above equation are plane waves: 

ip(x) = Ue ip ' x = B{p)e ip ' x , (5) 

where B{jp) is any function of the wave-vector p. 
Therefore, it follows from translational invariance 
that e lp ' x is an eigenfunction of □ with the cor¬ 
responding eigenvalue Bff) . Taking A to be a 
Lorentz transformation, it can be shown that Biff) 
can only depend on the the Lorentzian norm of p , 
i.e. p ■ p = rj^p^p 1 ', and whether or not p is future 
or past directed, i.e. sgn(p°): 

Bff) = B(sgn(p°),p ■ p ). (6) 

Combining Mj) and § we find B(—p) = B*(p ), 
which using Ml) is equivalent to 

Bffsgn(p°),p- p) = B{sgn(p°),p ■ p)*. (7) 

For a spacelike wave-vector p p , it is always possible 
to find a coordinate system in which p° = 0. As 
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a result, B(jp ) is real for spacelike p. For timelike 
momenta, however, B(p) may be complex and its 
imaginary part changes sign when p° —> —p°. 

Most nonlocal modifications of □ considered in the 
literature are of the form /(□), in which case B(p) 
is only a function of p -p. In this paper we focus on 
a class of nonlocal operators for which B(p) does 
depend on sgn(p°), and find many interesting con¬ 
sequences as a result. 

4. Locality at low energies: since □ provides a 
good description of nature at low energies, we re¬ 
quire □ —» □ in this regime. In other words, ex¬ 
panding B(sgn(p°),p • p) for “small” values of p-p : 
we require the leading order behaviour to be that 
of □: 

B(p)m%-p.p. ( 8 ) 

Note that by a “small” value of p ■ p, we mean in 
comparison to a scale which can be interpreted as 
the non-locality scale, implicitly defined through □. 

5. Stability: we require that evolution defined by □ 
is stable. This condition implies that B(p), when 
analytically continued to the complex plane of p , 
only has a zero at p ■ p = 0 [3] . 

6 . Retardedness: (D0)(j:) only depends on </>(?/), 
with y is in the causal past of x. 

Let us briefly consider a class of operators which satisfy 
all the aforementioned axioms. We shall let A denote the 
nonlocality energy scale and define 

A" 2 (n</>)(x) = acj){x) + A 4 [ f (A 2 t 2 )cf>(y)d 4 y, (9) 


Im(Z) 
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FIG. 1: Analytic structure of B(p) in the complex plane 
of Z = p ■ p/ A 2 



FIG. 2: The Fourier transform B(p) = g(p ■ p/ A 2 ) of □ 
defined in & where a and / are given by ©• 


where a is a dimensionless real number, J~ (x) denotes 
the causal past of X, and r xy is the Lorentzian distance 
between x and y: 

T xy = ( x ° — y 0 ) 2 — l x — y| 2 - (10) 

Examples of such operators have arisen in the causal set 
theory program EHS]. This operator is clearly linear, 
real, Poincare-invariant and retarded. It is shown in Ap¬ 
pendix [a] that there are choices of a and / for which □ 
is also stable and has the desired infrared behaviour ([ 8 ]). 
One such choice is 

4 e -s / 2 

f(s) = —5(s — e) --—(24 — 12s + s 2 ), a = —2, 

7T 47T 

(n) 

where e is an infinitesimally small positive number. 

The eigenvalues B(p) of □ take the form (see [3]) 


where p e is an infinitesimally small (p e ■ p e = — e 2 ), time- 
like, and future-directed (p° > 0) wave-vector. The an¬ 
alytic structure of B(p) is shown in Figure [I] Figure [2] 
shows the behaviour of B(p) as a function of p ■ p and 
sgn(p°) for the choice of / and a given in {mi¬ 


ni. CLASSICAL THEORY 

How would such non-local and retarded evolution man¬ 
ifest itself? To get a start on answering this question, we 
modify the evolution of a massless scalar field (j> coupled 
to a source J(x) via □ —» □: 

□0(x) = J(x) n<f>(x) = J(x). (14) 


A 2 B(p) = lim g((p + ip e ) • (p + ip e )/A 2 ), 
£-> 0 + 


( 12 ) 


/»00 

g(Z) = a + 4:nZ~^ / f(s 2 )s 2 K 1 {Z 1 / 2 s)ds, (13) 

J 0 


It is worth noting that the solutions of □</>(x) = 0 are 
identical to those of D(f>(x) = 0. This follows from requir¬ 
ing a stable evolution for □ (see [3]). As we will see in 
Section IIIB however, the story changes when J(x) 7 ^ 0. 
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A. Absence of an action principle 

It is natural to ask whether an action principle exists 
for (/>, whose variation would produce the non-local equa¬ 
tion of motion \3(f>(x) = J( x). One might propose to 
substitute □ with □ in the action of a massless scalar 
field: 


S[4>] = j d A x (^(f>{x)U<j>{x) - J(x)(f)(x) s j . (15) 

Requiring 5 , [0] to be stationary with respect to first order 
variations in <j> we find 3 

^(n + n T )<f>( x ) = j( x ), ( 19 ) 

where D T is defined in Fourier space via 

D T e ip ' x = B(p)*e ip ' x . (20) 


In the case of the retarded operator ([9]), for instance, 
0 T <j)(x) is the right hand side of with the domain 
of integration changed to the causal future of point x. 
Therefore, (15) does not lead to a retarded equation of 
motion. 

Due to the absence of a local Lagrangian description, 
quantizing a massless scalar field theory built upon □ is 


non-trivial. We shall address this problem in Section IV 


where we argue that the the Schwinger-Keldysh quan¬ 
tization scheme can still be used to obtain the desired 
non-local quantum field theory. 


B. Green’s function 

The Green’s functions of □ and □ are quite different, 
especially in the ultraviolet (UV) where their spectra dif¬ 
fer. One important difference is that □ , unlike □, has a 
unique inverse. Since □ is a retarded operator by defini¬ 
tion, it only has a retarded Green’s function. Recall that 


! To see this, it is instructive to express the action in Fourier 
space. Define the Fourier transform f(p) of f(x) via 

rf4p f(p)e ip x 


/(*) = / 


(16) 


Then, it can be shown that 


S = 


i 


d, A p 


</’(p)*^( s (p) + B{p)*)(t>(p) - (j>{p)*J{p) 


(17) 


Requiring S to be stationary with respect to first order variations 
4>{p) we find 


-(B(p) + -B(p)*)<Wp) = J{p)- 


(18) 


□ lias both a retarded G R (x,y) and advanced G A (x,y ) 
Green’s function: 

a x G R ’ A (x,y) = S^(x-y), (21) 

which satisfy the following “boundary conditions”: 
G R (x,y) vanishes unless x >- y (x is in the causal fu¬ 
ture of y), and G A (x,y ) vanishes unless y >- x. The 
two Green’s functions are related to one another via 
G A (x,y ) = G R (y, x). In the case of □, Green’s func¬ 
tion is unique (just the retarded one) and switching the 
arguments of the retarded Green’s function does not pro¬ 
duce another Green’s function. Let us show why this 
is. 

Let G(x,y) denote the Green’s function associated 
with □: 


U x G{x,y) = d (4) (x - y ), 


Note that G(x , y) can be expressed as 


G(x,y) 


d4 P 1 r iP-(x-v) 

(2tt) 4 B(p) 


( 22 ) 


(23) 


The path of integration in the complex p° plane is shown 
in Figure]^ This comes from the fact that □ is a retarded 
operator, so B(p) analytically continued to the complex 
p° plane takes its value above the cut. When B(jp) has 
no zeros in complex plane apart from at p ■ p = 0, which 
is guaranteed by the stability requirement, this choice 
of contour ensures that G(x,y ) = G R (x,y ) is indeed re¬ 
tarded. Switching the arguments of G R (x,y), we find 



d P 1 Jviv-x) 
(27r) 4 B[p) 

di P 1 r iy{x-y) 

(2tt) 4 B(-p) 

d P ^ Jp-G-v) 
(2tt) 4 B{p)* 


(24) 

(25) 

(26) 


where in the second line we have changed integration 
variables from p to —p. Then 



d 4 p B(p) 
(27t) 4 B(p)* 


pip-(x—y) ^ (5(4) (x - 


y)» 

(27) 


since B(p) is generically complex. As we will see in the 
sections to come, the fact that □ has a unique inverse 
plays a crucial role in the quantum theory of □. 


IV. QUANTUM THEORY 

We wish to construct a quantum theory of a mass¬ 
less scalar field <f> whose classical limit reproduces the 
retarded evolution induced by □. The quantization 
scheme which we believe is most suited in this case is the 
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FIG. 3: The integration path in the complex p° plane 
which defines the retarded G£een’s function associated 
with □. 


Schwinger-Keldysh (or double path integral) formalism. 
In what follows, we will first review the usual paths to 
quantization (i.e. canonical quantization and the Feyn¬ 
man path integral) and show why they fail in the case of a 
non-local and retarded operator like □. The goal of these 
discussions is to make clear why we choose the Schwinger- 
Keldysh formalism to construct a quantum field theory 
based on □. 


It is tempting to build the quantum theory of □ in a 
similar fashion: 


□</>(x) = 0 (32) 

[0(®)> Hv)\ = *A(x, y) = i(G R (x, y ) - G R (y, x)). (33) 

In this case, however, A (x,y) does not satisfy the equa¬ 
tion of motion (□ a; A(x, y) ^ 0) bec aus e G R (y, x) is not a 

and (|27|). There- 


Green’s function of □ (see Section 


fore, the equation of motion (321 is not consistent with 
the commutation relations (33). 


Ill 


It is worth noting that the root of this inconsistency is 
that the Fourier transform B(p) of □ is complex, which in 
turn follows from the fact that □ is retarded by definition. 
In Section IV C| we will arrive at a consistent quantum 
theory via the Schwinger-Keldysh formalism, using which 
we also build a Hilbert space representation of the theory. 
There we will see that the equation of motion (32) is 


given up in favour of the commutation relations (33). As 


it turns out, the degree to which (32) is violated depends 
on the imaginary part of B{jp). 


B. Feynman path integral 


A. Canonical quantization 

Let us consider the canonical quantization of a free 
massless scalar field <j>. The typical route to quantization 
is as follows: start from an action principal for cf> : derive 
the Hamiltonian in terms of 4> and its conjugate momen¬ 
tum, impose equal-time commutation relations, and fi¬ 
nally specify the dynamics via the Heisenberg equation. 
There is an equivalent approach, however, which defines 
the theory with no reference to an action principle, using 
the Klein-Gordon equation supplemented by the so-called 
Peierls form of the commutation relations: 

□^>( 2 :) = 0 (28) 

\cj){x),0{y)\=iA(x,y), (29) 

where A (x,y) is the Pauli-Jordan function: 


The Feynman path integral formalism requires a local 
Lagrangian description for the scalar field <j>. As was 
argued in Section [ill A[ however, this is not viable if one 
requires a retarded equation of motion. Therefore, the 
Feynman path integral formalism is also not suitable for 
quantizing this theory. 


C. Schwinger-Keldysh formalism 

The Schwinger-Keldysh formalism has a natural way 
of incorporating a retarded operator. In this approach an 
amplitude (called the decoherence functional T){(j ) + , <p~)) 
is assigned to a pair of paths {4> + ,<p~), which are con¬ 
strained to meet at the final time ( </> + (tf,x ) = (f>~ (tf, x)). 
The decoherence functional for a free massless scalar field 
takes the form 


A(x, y) = G r {x , y) - G A (x, y) 
= G R (x,y)-G R (y,x). 


(30) 


It is well known that (29) is entirely equivalent to, but 


more explicitly covariant than, the more commonly seen 
equal-time commutation relations (see e.g. Section C.2 
of [9]). Since A (x,y) is the difference of two Green’s 
functions, it satisfies the equation of motion: 


n x A(x,y) = 0. 


(31) 




where 


Exp 




d^x^ n R <j) cl + + ^(j) q n K <j) q 



(34) 

T (* + + r), 

(35) 


(36) 


This is why (28) and (291 are consistent with one another: 
both the left and right hand side of (29) vanish when 
is applied. 


In (34), is the retarded d’Alembertian, D" 4 = (D fl )^ 
is the advanced d’Alembertian, and D A is an anti- 
Hermitian operator which contains information about the 
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initial wave function m 4 5 Any source term J(a;) can be 
included by adding — J(f> + + J<j>~ = — \/2J(f> q to the inte¬ 
grand. 

Any n -point function in this theory is given by 

= J Dt+Dtjr^ixi) ■ ■ ■ ^ a ">(® n )z>(0+ 0"O37) 


where a, € {+,— ,q,cl}. These correlation functions are 
related to the correlation functions in Hilbert space rep¬ 
resentation by the following rule: 


( c/) + (an) • • • (f> + (x n )4> {yi) •••0 (y m )) 

= (0|T ^(yx) ■ ■ ■ T ${x{) ■ ■ ■ 4>(x n ) |0](38) 


where T{T) is the (anti) time-ordered operator, and |0) 
is the vacuum state of the free theory. 

In order to come up with a quantum theory for a non¬ 
local retarded operator, we replace D 77 with □ in (34) 
(and D K with \J K5 ). 


1. Classical limit 

Before going any further, let us take a look at the clas¬ 
sical limit of this theory. Performing Gaussian integrals 
(in the presence of a source term), we get 

{(f) c \x)) = -j=J d 4 y G R (x,y)J(y), (39) 
(<P q (x)) = 0, (40) 

resulting in 


-*< 

-*( 

-*< 


(j) cl (x)<t) q (y)) = G r {x , y) (42) 

(,b q (x)(f> cl (y )) = G a (x, y) = G R (y , x) (43) 

<j> c \x)<j> d {y))=G K {x,y) 


= - / d 4 zd 4 w G H (x,z)B K {z,w)G A (w,ylM) 


~i{ 4> q {x)<t> q {y)) = 0 


(45) 


where B K (x, y) is the kernel of □ /<6 . Using the definition 
of (f> q and (j ) cl , we get 


-*(</>+ {x)<j) + (y)) 


G K {x, y) + G R (x, y) + G A (x, y) , 

(46) 


-*(</> (x)<t> (y)) 


G K {x,y) - G R (x,y) - G A (x,y) , 

(47) 


0*0<a(2/)) = z G K (x,y)+ G H (x,y) - G A (x,y) 

( 4 “ 8 ) 

Note that if this theory has an equivalent representation 
in terms of field operator in a Hilbert space, then the 
above mentioned terms correspond to time-ordered two 
point function, anti time-ordered two point function and 
two point function respectively (see (38)). 

We require that the theory describes a free scalar field 
in flat space-time at its ground state. As a result, all n- 
point correlation functions of this theory must be trans¬ 
lation invariant, 


= (^ ai \ Xl +y) ■ ■ ■ ^\x n +y)). 

(49) 

This condition requires that all operators □, HT and \3 K 
must be translation invariant. Consequently, we get 


(0 + (z)) = (</> (x)) = J d 4 y G R (x,y)J{y). (41) 

It shows that in the classical limit where the field is rep¬ 
resented by its expectation value, there is no difference 
between </> + and (j>~ and both satisfy the retarded equa¬ 
tion of motion Ocf) = J. 


a K e i P -x = B K( p yi P x, ( 5 0 ) 

G K (x,y) = - j ^G R (p)B K (p)G A (p)e^-y\ 51) 

Note that D A is an anti-Hermitian operator. It 
means B K (p) is a total imaginary number (and 
G K (p ) = —G R (p)B K (p)G A (p) is also total imaginary 
since G R (p)G A (p) is real.) 


2. Green’s functions 

Let us consider the two point correlation functions of 
this theory in the absence of any source 


3. Fixing G K 

From here on, we assume that there is a Hilbert space 
representation of this theory with a Hamiltonian evolu¬ 
tion. We will justify this assumption later by finding 


4 The retarded and advanced d’Alembertians are defined via 
G R ’ A (D R ’ A f) = f for all suitable test functions /, where G R,A 
are the integral operators associated with the retarded and ad¬ 
vanced Green’s functions G R,A (x,y). 

5 We still need to determine \2 K . This has been done in 


IV C 3 


6 If & y (x) = <5! 4 )(x — y), then B K (x,y) = (Ci K S y )(x). With this 
definition, (□ K (j>){x) = f d 4 yB K (x, 
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the representation itself. In Appendix [B] we show that 
this assumption leads to the following relation, when the 
quantum system is in its ground state 


G K (p) = sgn(p") G H (p)~G A (p) 


(52) 


Note that (52) is nothing but the fluctuation dissipa¬ 
tion theorem (FDT) at zero temperature. This fixes the 
eigenvalues of D A as follows: 


B K (p) = 2ilml?(p)sgn(p u ) 


(53) 


4. Hilbert space representation 


We wish to find an equivalent Hilbert space represen¬ 
tation in terms of a field operator <j>(x) for this theory. 
As we mentioned earlier, (48) is the two point function 
of such a representation, 


W(x,y) = (0|^(x)^(y)|0) = (<j> (x)0 + (y)), (54) 


where |0) is the ground state. If we use (48) and (52), we 
arrive at 


W(x,y) 


d A p 2lm[B{p)\9(p°) ip . (x _ y) 

( 2 tt )4 \s(p)\2 ’ (55) 


where we call W(p) = 2 Im |^^| 2 ^ ^ ■ Since W(x,y) is a 

positive operator, Im[.B(p)]0(p o ) must be a non-negative 
number. So, we further assume 


sgn(Im[B(p)])=sgn(p°). (56) 

Once this condition is satisfied, the field operator <j>{x) 
and ground state | 0 ), defined to be 


m = /^\Aw(v ip ' x +sK ip ' x )-(57) 

[a p ,d q } = 5 {4 \p-q), (58) 

o P |0) = 0 Vp, (59) 

yield the desired correlation functions. 

Note that a p is only defined for time-like future- 
directed p, because otherwise W(p) is zero in the field 
expansion. It means that all tinre-like future-directed 
(positive energy) momenta contribute to the field expan¬ 
sion 157) . 


5. Hamiltonian 

By definition, time evolution operator is the operator 
that evolves <j>(x) in time, 

<J(t,x) = (60) 


It can be directly checked that 


U(t,to) = 

p-iHo(t-to) 

° ? 

(61) 

H 0 = 

J d A p p 0 d p a p , 

(62) 


gives the right time evolution. 

State |0) defined in (59) is the ground state of this 


Hamiltonian. Excited states (n-particle states) can be 
built by acting o^’s on | 0 ), 


\Pl ■■■Pn) = 


- a* 

X P1 


JO). 


(63) 


The excited state | p) represents a particle with energy p° 
and momentum p' where p° is independent of p 7 8 . This 
shows that the theory contains a continuum of massive 
particles with positive energy. The existence of a con¬ 
tinuum of massive particles in the context of Causal Set 
theory also has been pointed out in m, although their 
result is rather different in some other aspects. 


6. Comparison to local evolution 

At this point, it would be illustrative to consider the 
result of this formalism for LQFT. In this case 

B{p) = B local (p) = (p° + ie) 2 - |p| 2 , (64) 

where e is a small positive number taken to zero at the 
end of calculation. The two point function is given by 

w (p) = 2 (p 2)2 + (ep o )2 e (p°) = ^5(p 2 )0(p°). (65) 

As a result, 

W(x, y) = f jgfc2ir,6(p 2 )d(p°yr<*-y\ ( 66 ) 

£(*) = / 4k V 2 *Hp 2 Mp°) (a P e ip x + dj,e~ ip ' x ) (.67) 

Two point function and field expansion are exactly the 
ones we expected. Only on-shell particles {p ■ p = 0) 
contribute to the field expansion. 

Here, we see one important difference between local 
and retarded non-local evolution. In the local case, only 
on-shell modes (p ■ p = 0) contribute to the field expan¬ 
sion. As a result, excited states of the theory consist of 
all on-shell particles. In non-local retarded case (where 
generically Im[B(p)] ^ 0), off-shell modes (p-p/ 0) also 
contribute to the field expansion. Consequently, one ex¬ 
pect the existence of off-shell modes in ”in” and ’’out” 
state of scatterings in the interacting theory. 


7 Momentum operator P = f r/ 1 p p aj,a p is the generator of spacial 
translation. 

8 Note that these states are different from the usual states |p) used 
in LQFT which describe a particle with momentum p and energy 

IpI- 



















Let us investigate properties of W ( p ) for a generic non¬ 
local retarded operator. First of all, it is only non-zero for 
time-like future-directed momenta. This means that only 
time-like future-directed momenta contribute to the field 
expansion and can exist in ”in” and ’’out” state (particles 
with time-like momentum and positive energy). 

Considering that B(jp) is only zero at p -p = 0, W(p) is 
a finite number for all p-p 7^ 0 (we will see the significance 
of this result in IVI Bh . On the other hand, since in the 


subspace of on-shell modes □ operator is exactly the same 
as □, we conclude that W(p) = 2i tS(p 2 )9(p°) for p-p = 0. 
Therefore, W(p) consists of a divergent part at p ■ p = 0 
and a finite part for p ■ p ^ 0. This means that there are 
two different contributions to the field expansion (571, 


one from on-sliell modes that is the same as (67) and 
one from off-shell modes which only exists in the case of 
non-local retarded evolution 

£( X) = / (02 \/2 nS(jP)9(jP) (ape*-* + 


+ 


/p 2 /0 


d 4 p 

(2 W 2 


W{ P ) (a p e ip x +dle~* p - x ) . (68) 


then to define the two-point function to be the positive 
part of iA(x,y) in the following sense: 

(0|#z#(2/)|0) =Y / %T P (x)T;(y). (71) 


Taking G R (x , y) to be the retarded Green’s function of 
□ (see (23) and ( |26| ) , we find 


a a ^P‘ x 2Im(_£?(p)) ip. x 

,Ae =imr e ■ 


(72) 


which using the SJ formalism then leads to the two-point 
function 

(oi«x)?(»)io )=j 

(73) 

If condition (56) is satisfied, this two-point function is at 


that derived from the Schwinger-Keldysh formalism (see 
(55) and (56)). It is reassuring that two different paths 
to quantization, at least at the free level, lead to the same 
theory. 


D. Sorkin—Johnston quantization 

The Sorkin-Johnston (SJ) proposal defines a unique 
vacuum state for a free massive scalar field in an arbi¬ 
trarily curved spacetime m- This proposal is a contin¬ 
uum generalization of Johnston’s formulation of a free 
quantum scalar field theory on a background causal set 
m- As is the case for □, canonical quantization does not 
admit an obvious generalization for a causal set. The SJ 
quantization scheme uses only the retarded Green’s func¬ 
tion G R (x,y) to arrive at the quantum theory. Since □ 
also admits a retarded Green’s function, one can apply 
the SJ prescription to arrive at a free quantum field the¬ 
ory of the massless scalar field we have been considering. 
In what follows, we will show that the SJ proposal ap¬ 
plied to □ produces the same free quantum theory as the 
Schwinger Keldysh formalism, provided condition (56) is 
met. 

Consider the corresponding integral operator of the 
kernel iA(x,y) = G R (x,y) - G R (y,x): 


(iAf)(x)= / iA(x,y)f{y)d 4 y. 


(69) 


It can be shown that iA is Hermitian, which implies it 
has real eigenvalues, and that its non-zero eigenvalues 
come in positive and negative pairs: 


(iAT p )(x) = ApT p (x) 


(zAt;)(x) =-\ 2 p t;(x). 

m 

We have assumed here that the eigenfunctions T p form 
an orthonormal basis of L 2 , which can always be achieved 
since iA is Hermitian. The Sorkin-Johnston proposal is 


V. INTERACTING FIELD THEORY 

Let us now consider the interacting theory. We intro¬ 
duce the interaction in the Hilbert space representation 
by adding a potential term to the free Hamiltonian as 
follows: 


H(t) = H 0 + / d 3 x.V x)). 


(74) 


Starting with a general initial wave function, one is able 
to find the final state of the system by solving Heisenberg 
equation of motion in principle. However, in practice this 
is a very hard task to do. So, we try to find the S-matrix 
amplitudes perturbatively. 

In order to do so, we can use the available machinery 
of LQFT, and move to the interaction picture. Time 
evolution in the interaction picture is given by 


fj = j'ef di ' xV ($i) 


(75) 


where <pi is the held in the interaction picture given 


by (57). Perturbative expansion of t/j yields S-matrix 


amplitudes. Performing the calculations to find the S- 
matrix, we come up with modified Feynman rules for this 
theory. We explain these modifications in the following 
two examples. 


A. Example 1: 2-2 Scattering p\P 2 —> qiq 2 in 
theory 


Scattering amplitude S q 

'-’9192, P1P2 = (qiQ 2 \Tc 


is given by 




t \piP2)- (76) 
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To first order in A, it yields 


S qiq2 ,p lP2 = -i£fd 4 x (qiq 2 m(x)\PiP 2 ) 

= ]/w( Pl )W(p 2 )W(q 1 )W(q 2 )S( 4 HEP - E9), (77) 


where we have substituted for tfij from (57). It is inter¬ 
esting to note that ([77f is time reversal invariant. 

In the transition from local to retarded non-local prop¬ 
agation, here we see the first change in the scattering am¬ 
plitudes. The values assigned to each external line have 


changed from y/2n5(p 2 )9(p°) to y W(p). Note that here 
the scattering amplitude is computed in the basis of 4- 
momentum | p) which is different from 3-momentum basis 
|p) of LQFT. 


B. Example 2: 2-2 Scattering pip 2 —> qiq 2 in 2)</> 3 
theory 



FIG. 4: The amplitude of this diagram in LQFT is zero, 
because of the energy momentum conservation; two 
massless particles cannot produce a massless particle. 
However, in our theory there is a continuum of massive 
particles and the amplitude of this scattering is 
generically non-zero. 


In this case, S qiq2tPlP2 is given by 

S qiq2tPlP2 = (qiq 2 \Te-^ di ^\ PlP2 ). (78) 
To second order in A, it yields 
S qiq2lPlP2 = ^(=ir) 2 1 d 4 x d 4 y {q^Tcfc (a;)^(t/)|pip 2 ) 

= ^^W(p l )W(p 2 )W(q 1 )W(q 2 )5^(Y,P- E l) 

x G F (p\ + p 2 ) + G F (pi — q±) + G F (jp\ — q 2 ) 
G F (p) = i^y + 8 b~( p ) the time-ordered two point func¬ 


tion (46) in Fourier space. In the transition from local 


to non-local operator, here we see another change in the 
scattering amplitude. The values assigned to each inter¬ 
nal line have changed to the new value for the Feynman 
propagator Gf(p)- 

From these examples, it is obvious how scattering am¬ 
plitudes can be computed in this theory. For any Feyn¬ 
man diagram only the values assigned to external lines 
and internal lines have changed. Note that the amplitude 
of some diagrams in LQFT is zero, as a result of energy- 
momentum conservation, while in this theory they are 
not. For example in LQFT A cf> 3 theory, the amplitude 
assigned to diagram [4] is zero, because the sum of two 
(non-parallel) null vectors cannot be a null vector. How¬ 
ever, in this theory there is a continuum of massive par¬ 
ticles, and for example two on-shell particles can interact 
and produce one off-shell particle. 


then the differential transition rate of such scattering is 
zero. It means that in order to have a non-zero transi¬ 
tion rate (and cross-section), all of the incoming particles 
must be on-shell. This is the most distinctive property 
of off-shell particles: cross-section of any scattering with 
off-shell particles is zero. 

For now consider the scattering from state |a) = 
|pi • • • PN a ) to |/3) = \qi ■ ■ ■ qNp) where all the incoming 
particles are on-shell, pf = 0. Assuming that the inter¬ 
actions happen inside a box with volume V (see El), 
differential transition rate is given by 


dr = 2 tt 


AT 0 + l 


( 2 tt> 

V 


E, 


pi 


■E, 


Pn, 


x \Mp a \ 2 d 4 qi ■ ■ ■ d 4 qjyg, 


-5 (4) C^pi 

(79) 


where E Pi = |pi| and 


Sp a = -2m5^ 4) (y^ i p i -y^ j q i )\Jw(pi) • • -W(p Na )Mp a . 

(80) 

In the case of 2-2 scattering, the differential cross section 
is given by 

dcr = ^ d (4 \^p i -^q i )\ M p a \ 2 d 4 q 1 d 4 q 2l 

V X J Pi^P2 U 

(81) 

where 


VI. FROM SCATTERING AMPLITUDE TO 
TRANSITION RATE 


u = 


\/(Pl-P2) 2 -pIpI 


P1P 2 


(82) 


At this point, we want to find the rate of a process 
using the S-matrix amplitudes. In |VIB| we have shown 
that if one (or more) of the incoming particles is off-shell, 


is the speed of particle 1 in the frame of reference of 
particle 2 (and vice versa) and “ is the flux of incoming 
particles. 
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A. 


P1P2 


q±q 2 cross section in fi(j> 4 


where p 2 = jAp 2 . (871 is the same as (86) in the limit 
A/i 2 


0. 


As an example, we will find the cross section of P 1 P 2 —¥ 
qiq 2 where p 2 = 0. Using (77) and the definition (80), to 
first order in A 


The following field operator will yield the above two 
point function 


M = 


A 


W( qi )W(q 2 


(83) 


E 

j=i 


( 2 ? r) 5 

As a result, cross section is given by 

y _ _ 

da = — —-j-- -W(q 1 )W(q 2 )S { ‘ l \p 1 +p 2 -qi-q 2 )d 4 qid 4 q 2 . where 

4(2t r) \Pi-P2\ 

(84) 

Let us constraint the outgoing particles to be only on- 
shell q 2 = 0. In this case W functions in (841 pick up a 
delta function and one can check that ([84 ) for outgoing 
on-shell particles results in the usual cross section of X(f> 4 
in LQFT. However, if we constraint (at least) one of the 
outgoing particles to be off-shell with a fixed mass, the 
cross section becomes zero. Cross section over outgoing 
off-shell particles is only non-zero when the integration 
over continuum mass is also performed. We see the sig¬ 
nificance of this in the next section when considering the 
scattering of off-shell particles. Due to the contribution 
of off-shell states, the total cross section (84) is increased 
compared to the local theory. 


= / ( 2 % 

d 3 p 


ip-x 


+ c - c )Lo = 


p°=ipi 




(2^)3/2 x /2U^; 




P ip-x 


E p ,n = \/p 2 + P 2 


,pi 5 1 


Jt 

a 


= (5(3) (p - q )6 I 

p,m| 0)=0 


Pi, P-j 


+ C.c 


( 88 ) 

(89) 

(90) 


and state |p,/r) = |0) is a one particle state with 

momentum p, mass p and energy E p)|U . 

From now on, consider a concrete example of 2-2 scat¬ 
tering with Jy (jfi interaction and incoming particles with 
definite mass and momentum. The idea behind this proof 
can be generalized to more complicated examples. Up to 
first order in A 


Pi, tni; p 2 , TO 2 |<S'|qi, /xi; q 2 , p 2 ) = 
iX 


B. Off-shell particles and cross section 

In order to calculate the cross section of any scattering 
involving incoming off-shell particles, we make use of the 
fact that off-shell particles can be thought as a continuum 
of massive particles. 

This can be done by expressing the two-point function 
as a sum over massive two point functions: 


(2tt) : 


Zp-Z' 


\ * = 1 


AE E 

,pi -^Pi ,rrii 


In (91), if any of the particles was on-shell (say p\ =0), 
we should set IXp 2 p(p\) = 1, otherwise p is replaced by 

P- 

The differential cross section is given by 


W(x,y)= dp? p(p 2 ) 


d 4 p 

(2tt) 


da = (27r) 


-2 


4 2ird(p 0 )5(p-+p-)e 


2\ ip-(x-y) 


A 2 i A P' )' p(pi)p(pl)p(mi)p(m^) 

lQuEp 1 mi E p 2 _rjiy f 11 Uq 2 ^2 

5 (4) (pi +p 2 - qi - q 2 )d 3 pid 3 p 2 . (92) 


(85) 

where p(—p 2 ) = for p° > 0. Note from (68) that 
p(p 2 ) = S(p 2 ) + p(p 2 ) where p is a finite function. In 
other words, 


In order to get the total cross section, we should also sum 
over the mass parameter in the phase space of outgoing 
particles. In the (mass) continuum limit this means 


W(x 7 y) = 


d 4 p 


(2n) 4 


dp 2 p(p 2 ) 


2Tre(p°)6(p 2 )e ip{x - y) 
d 4 p 


E A ^ 2 p( m i) 


dm 2 p(m ■) 


(93) 


(2„)i 2M{ P 


2\ ip-(x-y 


P) 


} (. 86 ) 


In order to make everything more similar to LQFT, we 
discretize the mass parameter to get 


W(x,y) = 


d 4 p 


•EV^i) 

3 =1 


2ire(p°)5(p 2 y p d x - y) 

d 4 p 


(2ny 


-2n0(p")6(p 


p 2 j )e ip < x ~v\87) 


which absorbs two factor of A p 2 in (92); however, there 
are two remaining factors of A p 2 . If the incoming parti¬ 
cles (even one of them) are off-shell, since p(p 2 ) is a finite 
number, in the limit A p 2 —> 0, the cross section becomes 
zero. This means that the (total) transition rate of scat¬ 
tering with off-shell particles with fixed mass is zero. The 
cross section is only non-zero when both of the incoming 
particles are on-shell. 

This is, in fact, consistent with what we have found 
in the previous section. There, we have shown that the 
transition rate of on-shell —» off-shell is non-zero, only 
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when the integration over mass of the off-shell particles is 
performed. In fact, scattering transition rate of on-shell 
particles to off-shell particles with fixed masses is zero. 
Since the theory is time reversal invariant, this suggests 
that the scattering transition rate of off-shell particles 
with fixed masses must be zero too; consistent with what 
we have found here. 

This also means that an initial state with a suitable 
continuum superposition of off-shell masses can scatter 
into on-shell modes (time reverse of the process of on- 
shell scattering into off-shell). However, as we argue in 
the next Section, these states are fine-tuned and generally 
we do not expect to find the system in these superposi¬ 
tions. 


C. Off-shell —> on-shell scattering: continued 


In the previous section, we showed that the transition 
rate of scattering with off-shell particle(s) is zero. How¬ 
ever, a suitable continuum superposition of off-shell par¬ 
ticles can scatter non-trivially. In this section, we want 
to explain this point to a greater extent and argue that 
it is unlikely to find the system in these superpositions. 
We will not go through the detail of calculations since it 
is not essential to our argument in this secion. 

We make use of the following toy model theory that 
mimics many properties of the proposed nonlocal theory: 

1 N 1 

C = -V’onV’o + 2- AV' 4 , (9 4 ) 

i=l Z 
N 


This is a theory of one massless scalar field (playing the 
role of on-shell modes) in addition to N massive scalar 
fields (playing the role of off-shell modes) and we are in¬ 
terested in N —> oo limit of the theory (A and gi’s are cou¬ 
pling constants and do not scale with N). The advantage 
of working with this theory is that while its behaviour is 
very similar to the non-local theory, (94) is a local quan¬ 
tum field theory and possibly more comprehensible to 
the reader. The interaction term in (94) is designed in 
a way that interactions with massive (off-sliell) fields are 


suppressed by a factor of \/N and in N —> oo limit their 
interactions become negligible. On the other hand, the 
number of off-shell fields goes to infinity. In what follows, 
we explain that this theory imitate many properties of 
off-shell and on-shell particles in the non-local theory. 
First, let us define the following quantities: 

the scattering cross section of two 
particles with masses and momenta mi,pi and m 2 ,pi 
into two particles with masses pi and ^2 (VVi an d Vv* 2 ) 
and cT^fm 2 i s the total scattering cross section of two 
particles with masses and momenta mi,pi and 7772, j?2- 
Consider the scattering of two t/>o particles into two 
final particles. If we restrict the two final particles to 


be massive (off-shell fields with fixed masses), then the 
scattering cross section in N —> 00 limit goes to zero. 
However, if we sum over all massive final states (all off- 
shell particles), the total cross section is non-zero. In fact, 
for different final states the corresponding cross sections 
scales with N as follows: 


P 1 P 2 

^oo-^oo 

oc 

N°, 

-P 1 P 2 

oc 

1 

° 00—>-0m 

N’ 

.P 1 P 2 

oc 

1 

00— >mirri2 

N*’ 


m/0, 


777-1, m 2 yf 0. 


While the interactions with individual massive fields are 
suppressed, the number of massive states scales with N. 
In this way, the total scattering cross section of two initial 
massless particles into two massive final states, summed 
over all masses, is finite and non-zero (the same scaling 
works for scattering into one massless and one massive 
particles). 

On the other hand, any scattering with (at least) 
one massive initial state result into zero cross section. 
For example, the following total scattering cross sections 
(summed over all final states) scale with N as 


-P1P2 

°0m 

0 " 

db 

\k 

8 

(95) 

-P 1 P 2 
rap 7712 

OC -^ 2 , 7771,7772 ^0, 

(96) 


and they vanish in N —► 00 limit. 


As we showed, massive particles in this theory (94) 


mimic the properties of off-shell states in the non-local 
theory; they can be produced by the scattering of mass¬ 
less states, while the reverse process (scattering of mas¬ 
sive states into massless) does not happen. 

However, the theory is (obviously) time reversal invari¬ 
ant and massive massless scatterings must take place. 
This is indeed true, but as we demonstrate here the ini¬ 
tial massive state that scatters non-trivially must be a 
superposition of different masses. Consider state 7 , a su¬ 
perposition of M different masses, scatters off a massless 
particle. Then, the total transition probability T 07 scales 


as 


<A n 


(97) 


Where A have no dependence on M and N (see Appendix 
[C] for proof). This transition probability is non-zero in 
N —► 00 limit, only when M also scales with N. 

So, massive —> massless scattering indeed happens. 
However, the massive state that scatters non-trivially 
must be a superposition of (infinitely) many different 
masses and in this sense is fine-tuned. It is similar to an 
egg that smashes into pieces upon falling on the ground; 
the reverse process of pieces assembling an egg can in 
principle happen, but it is very unlikely. 

In this sense, we expect the off-shell to on-shell scat¬ 
tering in the non-local theory to be negligible. In prin¬ 
ciple this transition can happen, but it is very implausi¬ 
ble. The essence of our reasoning in this section is based 
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on thermodynamical arguments and although it is not a 
complete proof, we hope that we have provided enough 
evidence to show that off-shell —> on-shell scattering is 
very unlikely. Definitely, further quantitative studies are 
needed to augment (or disprove) our claim. Perhaps, a 
good starting point is to consider the toy model theory 
(94), since it shares a lot of properties of the non-local 
theory. 


Conditions 4 and 5 in Section |TT] and (56) must be re¬ 
placed by the above-mentioned conditions. One easy way 
to come up with such an operator is to make use of the 
existing operator B(p) in the massless case, and consider 
it as a function of p 2 and sgn(p°). Then, 


B m {p) = B(p 2 + to 2 , sgn{p°)). (99) 

has all the desired properties (this also has been shown 

in EU). 


VII. EXTENSION TO MASSIVE SCALAR 
FIELDS 


VIII. CONCLUSION 


Throughout the paper, we only considered the modifi¬ 
cation of a massless scalar field. But what about massive 
scalar fields? One may suggest to replace □ with □ in 
the equation of motion of a massive scalar field as follows 

(□ - M 2 )(j>(x) = J(x) (98) 


and follow similar steps of quantization. However, this 
method does not work. If M is a real number, then there 
is no mode satisfying ( [98]) in the absence of J. In other 
words, there is no on-shell modes. 

Another way is to choose M to be a complex number 
such that for a time-like future directed momentum p 1 
B{p) = M 2 . In this case, the mass of on-shell mode is 
given by m 2 = — p 2 . However, □ — M 2 is no longer a real 
operator and the solution to (98) generically cannot be 
real. 

The extension to massive scalar fields can be done by 
considering the following observation. All of the proper¬ 
ties in massless case can be read from the analytic struc¬ 
ture of B(p) in Figure [TJ Massless modes are on-shell 
because there is a simple zero at p 2 = 0 and there are 
off-shell modes for time-like momenta because there is a 
cut for time-like momenta in [T] 

In this way, the extension to massive case seems much 
simpler. □ — m 2 must be replaced with O m whose eigen¬ 
values B m (p ) satisfy the followings: 


1. There is only one simple zero at p 2 = —to 2 . Also 
lim p 2 +rra 2_ >0 = —1 to get the correct local 

limit. 


2. The cut must be only on momenta with higher 
masses p 2 < —to 2 . Otherwise, in the quantum 
theory, there are off-shell modes with mass smaller 
than to which makes the on-shell mode unsta¬ 
ble (on-shell modes can always decay into off-shell 
modes with less mass). 

3. Im B m (p) > 0 for p° > 0. 


In this paper, we studied the physical consequences 
of a causal non-local evolution of a massless scalar field. 
We started by modifying the d’Alembertian to a causal 
non-local operator at high energies. Quantization of a 
free field showed that the field represents a continuum of 
massive particles. In fact, there were two sets of modes: 
on-shell modes (massless particles) and off-shell modes 
(massive particles). 

The Feynman rules for the perturbative calculation of 
S-matrix amplitudes were discussed. The most impor¬ 
tant result (in our opinion) is the fact that the cross sec¬ 
tion of any scattering with off-shell particles is zero. This 
suggests that although these modes exist and probably 
can be detected by other means, there is no way of detect¬ 
ing them through scattering experiments. This property 
opens up the possibility that dark matter particles might 
be just the off-shell modes of known matter. Finally, we 
extended this formalism to massive scalar fields. 

Throughout this paper we only considered scalar field 
theories, but how about other types of fields? Exten¬ 
sion to other types of fields, such as vector held, is not 
as straightforward as for scalar fields. Incorporating 
gauge symmetry in the theory is another important issue. 
Whether causal Lorentzian evolution can be extended to 
vector fields (and other types fields rather than scalars) 
can be the subject of future studies. 
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in m • It follows from their result that 

g(z) = g(z-z), (A6) 

/»00 

g(Z) =a + 47rZ-5 / f(s 2 )s 2 K 1 (Z 1 / 2 s)ds , (A7) 

Jo 

where an infinitesimal time-like and future-directed imag- 
ina ry p art ought to be added to z on the right hand side 
of (|A6|) (see [3] for more details). 


1. IR conditions 


The infrared condition 


g(z) ^ -z. 


is equivalent to satisfying 

(AS) 


In [3], a framework is developed to determine what con¬ 
straints (A8) places on a and /, for some specific choices 


of / which arise in causal set theory. Generalizing that 
methodology in a straightforward manner, we find that 
(A8) is true if and only if the following conditions are 


satisfied: 


f(s 2 )s 2k+1 ds = 0, 


pOO 

/ /(s 2 )s 5 In sds = — 

Jo 

poo 

a + 27r / f(s 2 )s 3 In sds = 0. 

Jo 


k = 0,1,2 (A9) 

(A10) 
(All) 


Appendix A: Existence and Examples of □ 

Here we will show there are operators □ which satisfy 
all the axioms introduced in Section [IlJ In fact, we will 
outline a procedure for constructing such operators. 

We shall consider the following operator: 

A ~ 2 (□<?!>) (a;) = acfix) + A 4 f f{A 2 T 2 )(j){y)d 4 y, 

Jj-Q) 

(Al) 

where A denotes the nonlocality energy scale, a is a 
dimension-less real number, J~{x) denotes the causal 
past of x, and r xy is the Lorentzian distance between 
x and y: 



T xy = ( x ° ~ V 0 ) 2 ~ l x ~ y| 2 - 

(A2) 

It may be 

shown that 


Ue ip ' x 

= B(p)e 1 ^, 

(A3) 

B (p) 

= A 2 g(p/A), 

(A4) 

g(z) 

= a+ [ f((y°) 2 -\y\ 2 )e- iz - y d 4 y, 

(A5) 


Jj+( 0 ) 


where as 

usual x ■ y = rj^x^y 11 . Evaluating 

g(z) 


amounts to computing the Laplace transform of a re¬ 
tarded, Lorentz invariant function, which has been done 


2. From B(p) to □ 


It is often desirable to constrain the behaviour of B{p), 
as opposed to □ directly. For instance, as is argued in 
Section |IV C 4[ the quantum theory is well behaved only 
when the imaginary part of B(p ) (for timelike and future- 
directed p) is always positive. The question then be¬ 
comes: are there any choices of a and f which allow for 
this possibility, provided the IR conditions (A9) (All) 


are satisfied ? To answer this question, we turn the prob¬ 
lem around. Given a choice of B(p), we reconstruct a 
and / and then ask if the IR conditions are met. 

It can be shown that for x > 0: (see e.g. 10.27.9 and 
10.27.10 of US]) 

g{-x 2 - it) = gni-x 2 - ie) + igi{-x 2 - ie), (A12) 
Ojr r°° 

gB,{—x 2 — ie)=a-\ -/ f(s 2 )s 2 Yi{xs)ds, (A13) 

x Jo 

2tt 2 r °° 

gi(-x 2 —ie) = —-— / f(s 2 )s 2 Ji{xs)ds. (A14) 
x Jo 

We can now use the following orthonormality conditions 
of Bessel functions (see e.g. 1.17.13 of 023 ) to express / 
in terms of gj: 

poo 

5{x — x)=x / tJi(xt)J\(xt)dt. (A15) 

Jo 
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Doing so yields: 

f(s 2 ) = f g (s 2 )+h(s 2 ), 


f a (s 2 ) = 


1 


gi{—x 2 — ie)x 2 Ji(sx)dx, (A17) 


2 7T 2 S 

where h satisfies for all x: 

pOO 

/ h(s 2 )s 2 Ji(xs)ds = 0. 


The IR conditions then reduce to 

r°° /*°< 

( Al6 ) / 5/ (-z 2 - ie)x dx = 7T 2 / 

JO JO 

/*°° _ 7T 

/ 9i(-x 2 - ie)x~ 3 dx = -, 

/ gi(—x 2 —ie)x~ 1 dx = ——a. 

Jo 2 


h(u)du, 


(A 18) 


(A30) 

(A31) 

(A32) 


This means that specifying gj(—x 2 — ie) fix es / up to any 
part for which the right hand side of (A14) vanishes. One 
example of a nontrivial function which satisfies (A18| is 
the delta function: h(x) = 5 + (s) = S(x — e), where e is 
an arbitrarily small positive real number. 

We can now express the IR conditions in terms of gi 
and h: 


» 2 ^ 2fc+1 d s 


/ h(s 2 )s 

Jo 

| pOO pOO 

— —— J gi(—x 2 — ie)x 2 J dss 2k J\{xs) = 0 

pOO 

/ h(s 2 )s 5 In sds 

Jo using (|A30|)-(|A32|): 

r°° r°° 4 

/ gj(—x 2 —ie)x 2 / dss 4 Ji(xs) In s = -,(A20) 

Jo Jo n 


Note that the only nontrivial condition to satisfy is 
(A30), since (A311 just fixes the normalization of gi and 
(A32) determines a. Note that for positive gi(—x 2 — ie) 
which is required by consistent quantum theory, a must 
be a negative number. 

If h is taken to be zero, then gi ought to change sign, 
which leads to a quantum theory with an unbounded 
Hamiltonian. We note that the class of operators which 
arise in causal set theory in [3] all have h = 0, and there¬ 
fore this feature. 

Let us work out a complete example in 4D. Let 
(A19) gi(_x — ie) = Ax 2 e~ x ^ 2 , h(x) = a5 + (x). (A33) 

where A and a are real constants. It can then be shown 


1 

2-7T 2 


A = — 
2 ’ 


4 

a = 

7T 


a = —2. 


+27T / h(s 2 )s 3 hi sds 


1 
7r 


gi(-x 2 - ie )a 


dss 2 Ji(xs) Ins = 0. 


It then follows from (A17|) that 
(A21) fg(s) = ~- 


o-s/2 


The above integrals over s are not absolutely convergent, 
so use the usual trick: 


Therefore: 


) -| 

dsJi(xs)e~ Ss —°> —, 
x 


/»oo 

/ dss 2 Ji(xs)e 
jo 

/»oo 

/ dss 4 J\(xs)e 


— <5s >0 


35 

> 7' 


-Ss s^o k -455 


/»00 

/ dss 2 Ji(xs) lnse -,5s — 2x~ 3 , 

Jo 


dss 4 J±(xs )In se 


—Ss ^0 


> 16x 


-5 


(A22) 

(A23) 

(A24) 

(A25) 

(A26) 


47T 

4 , e -s/2 

/(») = -> (-) - -Jf 


(24- 12s + s 2 ). 


-(24- 12s + s 2 ). 


(A34) 


(A35) 


(A36) 


Having the delta function example in mind, we shall re¬ 
quire h to satisfy for all k = 1, 2 


h(s 2 )i 


/»00 

l ds = 0, / h(s 2 )s 2k+1 In sds = 0, 

Jo 


3. Stability from positivity of gi 

We have required that evolution defined by □ should 
be stable. Instabilities are in general associated with “un¬ 
stable modes”, and in line with [3], we shall use this as our 
criterion of instability. More specifically, we take such a 
mode to be a plane-wave e lp ' x satisfying the equation of 
motion De lp ' x = 0, with the wave-vector p possessing a 
future-directed timelike imaginary part (i.e. p = PR + ipi 
where pi - pi < 0 and p 3 > 0). It is shown in [3] that the 
necessary and sufficient condition for avoiding unstable 
inodes is 


g(Z) ^ 0 , V Z ± 0 and Z e 


(A37) 


(A27) 

Also, we assume that the following integrals converge: 


pOO 

/ gi(—x 2 — ie)x~ k dx 

Jo 

pOO 

/ gi{—x 2 — ie)x~ k \nxdx 

Jo 


< 00 , 


< oo 


k = 1,2,3,4 

(A28) 

A; = 2,4. (A29) 


On the other hand, we argued in |IV C 4] that for consis¬ 
tency reasons we need to assume ImB(p) > 0 for p° > 0 
which implies g(Z) has a positive (negative) imaginary 
part under (above) the cut in Figure [l] 

Here, we show that not even stabilit y co ndition and 
positivity of gi(—x 2 — ie) (see Appendix | A 2[ ) are consis¬ 
tent, but latter is a sufficient condition for stability. In 
order to prove it, we make the following assumptions: 
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FIG. 5: The integration path in the complex Z plane. 
The closed contour is taken to be counterclockwise. 


1. g(Z) has a simple zero at Z = 0. IR conditions on 
g{Z) ( |A 8 [ ) guarantee this assumption. 

2. g(Z) has positive (negative) imaginary part under 
(above) the cut. 

We prove this by counting the number of zeros of g 
inside contour C = C\ + C 2 + C 3 + C 4 in Figure [5] 

If N and P are the number of zeros and poles of g , 
respectively, inside the contour C (taken to be anticlock¬ 
wise), then 


9'{Z) 


lc dZ g{z) 


= -2? ri(N-P). 


(A38) 


Let’s evaluate the left hand side of (A38 1 for each contour 
separately: 

1. According to ( |A7| , g(Z) approaches the con¬ 
stant value of a < 0 (see A2) for large Z. In fact, 
g(Z) —>• a + G(jpr) for some positive value of n 
(which depends on the function /). This means for 
« / 0 , 


X/ Z ii= 0 ' (A39) 

2. C 2 & C 4 : Since the values of g above and under the 
cut are complex conjugate of each other, the contri¬ 
bution from these diagrams can be added together 
to get 


9'{Z) 


/ dZ 

IC 2 +C 4 9\Z) 


Mm I" d X 9 ' {x + ie) 


g(x + ie) 


= 2 ilm In 


9(0 + ie) 
g{~ 00 + ie) 


,(A40) 


where e is an inhnitesimal positive number. 

If we define g(Z) = r g (Z)e l(p »^ z \ the right hand 
side of ( A40| ) (apart from the factor of 2 i) measures 
how much ip g rotates from Z = —00 + ie to Z = 


0 T ie. Since Img(x + ie) < 0 on the whole negative 
real line, In [g(x + ie)] is definable on one Riemann 
sheet. Combining this result with g {—00 + ie) = 
a < 0 and g {0 + ie) = — ie , we get 


f ■ 

/ dZ ~TZ\ = l7T - 
IC 2 +C 4 9\Z) 


(A41) 


3. C 3 : IR conditions require that close to Z = 0, 
g{Z) = —Z. This means 



9\Z) 

9(Z) 



1 

Z 


= — in. 


(A42) 


Adding the values of all the contours and considering the 
fact that g(Z) is finite everywhere {P = 0), we conclude 
that the number of zeros of g in complex plane of Z 
(inside contour C) is zero. Since there is no zero on the 
negative real line (Img(.T + ie) 7 ^ 0), there is no zero of 
g in the complex plane of Z except the one at Z = 0. 
Therefore, stability has been proven. 


Appendix B: FDT 


Here, we present the proof of (52) 9 . Let’s start by the 
following definitions 


iA(x,y) 

III 

' ^>' 
8 ^ 

5=P 

(Bl) 

G (1 \x,y) 

= ^{<£(x),<£(y)j^ , 

(B2) 

W + (x,y) 

III 

xix 

(B3) 

W~ (x,y) 

= ($(y)Hx )), 

(B4) 

G F (x,y) 

53 

(h) 

h 

•<?* 

1 

III 

(B5) 

is anti-commutator and () shows 
quantum state. If we define 

expectation 

G R {x,y) 

= A(x,y)H{x y y), 

(B 6 ) 

G A {x,y) = 

-A(x,y)H(x -< y ), 10 

(B7) 


we get the following relations 

iA(x,y) = W + (x,y) — W~(x,y) 

= i [G fl (z, y) - G A (x, y)] , (B 8 ) 

G {1 \x,y) = W + (x,y) + W~(x,y), (B9) 

G A (x,y) =G R (y,x), (B10) 

G F (x, y) = ^[G R {x,y) + G A (x,y)\ - ^G {1) (x,y)( Bll) 


9 Most of the content of this appendix is taken from m- 
10 where H is the Heaviside function: H(x y y) = 1 if x y y and 
otherwise 0 
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For a translational invariant system, the value of all the 
two point functions depend only on space-time separa¬ 
tion. This will allow us to define the following Fourier 
transform with respect to time 

A(cj,x,x / )= f dt A(t,x\t',x!)e~ lul ^ t ~ t \ (B12) 


Now, let us assume that the quantum system is in ther¬ 
mal state with temperature T = j ,. It requires that 

W ± (t, x; t ’, x') = W T (t + i/3, x; t' (B13) 


resulting in the following relation in Fourier space 

W + (co,x,y)=e^W~(u J ,x,y). (B14) 


Using (B 8 ), we get 


W%,x, y) 


W (w,x,y) 


iA(u>, x, y) 

1 - ’ 
x, y) 

1 - 


(B15) 

(B16) 


which reduces to 

Append 


(52) at zero temperature. 

ix C: Quantum Transition 


We start by proving a simple theorem for any quan¬ 
tum system. Consider a quantum mechanical system in 
the (normalized) initial state |a) evolves in time and the 
probability of finding the system at a later time tf in the 
state | fii) is called Pi, and assume |/3j)’s are orthonormal: 

Pi = \(/3i\U\a)\ 2 (Cl) 

where U is the time evolution operator. 

Now, consider a (normalized) state |/3) as a superposi¬ 
tion of \/3i) states: 


|/?> = E C ‘IA> ( C2 ) 

i 

£n 2 = l 

i 


Probability P of measuring the system at time tf in the 
state |/3) is given by 


On the other hand, since G R and G A are time transpose 
of each other, in Fourier space they are complex conju¬ 
gate. As a result, 

ImG F (w,x,y) = -^ReG ( 1 ) (w,x,y) 

= -1 [W + (co,x,y) + W~(u;,x,y) 

= -^*A(w,x,y)coth(^) (B17) 


where Im and Re are imaginary part and real part respec¬ 
tively and in the second line we have used the positivity 
of two point function W + (resulting that W (w,x, y) 
and W (w,x, y) are real.) 

With the assumption that this field theory in Hilbert 
space representation has an equivalent representation in 
terms of double path integral, time ordered two point 
function is given by (|46[). In Fourier space, it reads 


d (w, x, y) = - \cf (w, x, y) + (f (w, x, y) + (f (w, x, y) 

(B18) 

(d (w, x, y) is a total imaginary number and 

G^uqx, y) + G^oqx, y) is a real number. As a 
result, 


G* (w, x, y) = 2iha(f (w, x, y). (B19) 


Then, 


P = Mu\u)\ 2 - 


(C3) 


P =|(/3|[/|a )| 2 = 


E c *wi a > 

i 

< (em 2 ) (Ek&IWI 2 ) 


= E* 


(C4) 


where we have used the triangular inequality in the sec¬ 
ond line. So P is bounded from above by JV Pi. 


Now, let’s get back to the scattering of a massless par¬ 
ticle with state 7 , a superposition of M different masses, 
in Section VIC We already have shown (see (95)) that 
romi defined as transition probability of a massless par¬ 
ticle scattering with a massive particle (mass ?rq) scales 
with N as 


p _ _i 

1 0m ' ~ N 


(C5) 


where Ai depends on the momentum of the particles but 
independent of N. Using (C4| for transition probabili¬ 
ties, we conclude that 


Combining ( |B 8 | |Bl7l ) ( [bT 9] ) we arrive at 
G K (u,x,y) = coth( ~) [g^w, x, y) - (f (u, x, y) 


v—' Ai Ad 

i 

(B20) where A is the maximum of A/s. 


(C 6 ) 





















